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Abstract 

The visibility times and communication link dynamics for the EOS-McMurdo Sound direct downlinks have 
been studied recently. The 16-day EOS periodicity may be shown with the Goddard Trajectory 
Determination System (GTDS), and the entire 16-day period should be simulated for representative link 
statistics. We desire many attributes of the downlink, however, and a faster orbital determination method is 
desirable. We use the method of osculating elements for speed and accuracy in simulating the EOS orbit 
The accuracy of the method of osculating elements is demonstrated by closely reproducing the observed 
16-day Landsat periodicity. An autocorrelation function method is used to show the correlation spike at 
16 days. 

The entire 16-day record of passes over McMurdo Sound is then used to generate statistics for innage time, 
outage time, elevation angle, antenna angle rates, and propagation loss. The elevation angle probability 
density function is compared with a 1967 analytic approximation which has been used for medium to high 
altitude satellites. One practical result of this comparison is seen to be the rare occurrence of zenith passes. 
The new result is functionally different than the earlier result, with a heavy emphasis on low elevation 
angles. EOS is one of a large class of sun synchronous satellites which may be downlinked to McMurdo 
Sound. We examine delay statistics for an entire group of sun synchronous satellites ranging from 400 km to 
1000 km altitude. Outage probability density function results are presented three dimensionally. 

1. Introduction 

A new series of environmental spacecraft, such as the Earth Observation Satellite (EOS), may use the Tracking 
and Data Relay Satellite System (TDRSS) in some unique ways. Polar orbiting environmental spacecraft spend an 
appreciable portion of their orbit over the poles. The polar orbit of EOS or of Landsat may be represented as in 
Figure 1-1. The individual views are taken as a series of snapshots at 8.2 minute intervals. It is seen to be slightly 
retrograde, with a 98-degree inclination. EOS and other environmental spacecraft may find it more convenient to 
downlink their data directly to a ground station in Antarctica rather than relaying it over a distance of six earth 
radii to a geosynchronous relay. We examine these downlinks in some detail. 

A ground station at McMurdo Sound, Antarctica has been constructed by NASA to answer these needs. A polar 
satellite such as an Earth Observation Satellite (EOS) would be able to link to McMurdo Sound in a visibility 
region as shown in Figure 1-2. McMurdo would then uplink to TDRS West, to complete the link to the White 
Sands Complex (WSC). 
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Figure 1-1: Inclined Landsat (EOS) Orbit Inside VanAllen Belt 



Note the fint EOS p*u, labelled (1). 

Figure 1 - 2 : McMurdo Sound-EOS Visibility Region Rm/iut 25.6° 
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We address some of the communication link dynamics for the user satellite pass over McMurdo Sound. Sun 
synchronous satellites will be of widespread interest, and we closely examine arrival times, departure times, and 
atmospheric signal loss as a function of time during the pass over the ground station. These data may be further 
examined to generate interarrival time statistics, link duration statistics, and signal loss statistics as an aid to the 
system designer. EOS, as a clearly important member of the class of sun synchronous satellites, gets special 
attention. 

Many long computer simulations are necessary for this report. The method of osculating orbits (1) is used here for 
speed and accuracy. The oblate earth perturbation terms are vital to descriptions of the 16-day EOS periodicity. 
Indeed, the 16-day periodicity can be used as a simulation check against existing Landsat4 data. Most of our EOS 
runs include the entire 16-day periodic interval as a key check. However, only small portions of the output tables 
are included here because of space limits. The report begins with a brief description of oblate earth perturbation 
terms which are necessary for the method of osculating elements. 

Atmospheric attenuation on the McMurdo-satellite link is closely related to the orbital computations because the 
attenuation depends strongly on the ground station elevation angle. We use two new attenuation results as an 
integral part of the orbital program. The fust new propagation result is due to Liebe (2). The new propagation 
results are used but not included because of space limitations. A short portion of a long EOS simulation follows. 
The long-term EOS statistics are included as a way to keep the results brief. General sun synchronous results 
conclude the report. Sun synchronous altitudes between 400-1000 km will be seen to offer short interarrival times 
at McMurdo Sound, but altitudes below 500 km will occasionally skip an orbit before visibility repeats. 


Appendix A condenses the long 16 day simulation into short graphical output. 

2.0 Analysis 

The anticipated long 16-day periodicity of EOS leads us to consider methods to reproduce the entire 16-day period 
with suitable accuracy. Anything less would cause serious doubts about the ground elevation statistics and the 
propagation loss studies which follow. The method of osculating elements, with analytic approximations for the 
perturbed motion, is useful for fast and relatively accurate orbit prediction for missions of long duration. Moulton 
(3) gives the basis of the method as the partial derivatives of the instantaneous orbital elements as a function of a 
perturbing potential. Brouwer (4), Kozai (5), and Blitzer (6) have used the method to achieve exhaustive insights 
into the motion of low earth orbiting satellites. Their results for satellite perturbations may be subdivided into 
secular components, long cycle components, and short cycle components. We desire only the secular components 
here, and omit the long and short cycle terms because they tend to cancel out their effects over the long term. 
Omission of the short cycle component has implied a maximum instantaneous orbital error of a few kilometers for 
the 1958 Vanguard orbit. We accept that error, along with the error due to lunar perturbations (although lunar 
perturbations can be included with the method of osculating elements with the aid of results by Ash (1)). 


The key secular terms are: 




• ^4 + 12JT-7 1 - 9e 3 )cos(i) + (-40- 36-J7--7 2 +5e’ )cos’(i)] 
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with the key variables defined below 

The method of osculating orbital elements allows the Keplerian orbital elements to be continuously updated 
without concern for the time interval between updates. It also offers a concise orbital description, which allows 
other important satellite communication features such as link propagation to be included in the same computer 
code. 

The secular perturbations may be abbreviated further, while usually retaining over 99% of the accuracy of the 
secular terms. The oblateness effects have been found to give a regression of nodes 


W =- 


2a 



cos( i) rad/se c 


2-4 


The argument of perigee changes as 


W r 



(-1 + 5cos 2 (i)) rad/ sec 


and the mean anomaly changes as 



<J}-e 2 (-/ + 3cos 2 (i)) radj sec 


2-5 


2-6 


W = right ascension, rad 
Wp = argument of perigee, rad 
R = earth radius, km 
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Where 


p = earth gravitational constant 
= 0.39860064 * 10 6 km 3 /sec 2 

P = a (l-e 2 ) = semilatus rectum 

h = 1.082635 * 10' 3 

e = eccentricity 

i = inclination, rad. 


3.0 EOS Simulation Results 

Table 3-1 indicates a first pass (labeled 1) of EOS over the McMurdo Sound ground station. It uses 1992 Landsat4 
orbital elements. Table 3-1 lists the first and second passes in an abbreviated simulation output. Elevation angle, 
attenuation, and subsatellite point are listed as a function of time for the EOS passes. The first attenuation column 
is for a single link, and the second represents a dual switched link. Attenuation is chosen at the 0.997 precipitation 
probability level for a 1.0 km altitude precipitation layer (7,9). The subsatellite point is recognized as North 
Latitude and East Longitude. Time is in hours from start of the Epoch (note 1 minute increments). Range is in 
megameters. 


Table 3-1: EOS Downlinked to Antarctica 

EOS DOWNLINKED TO ANTARCTICA TAVG(K) * 201 

EOS LINK ; X BAND GND SITE-78.11 166.2272 GFHZ* 8.2875E+09 

RP PREC. RATE (MM/HR) * 1.4 98039 AVAIL* .997 PRECIP. HT.= 1 
A,E,I* 7077.79 .0014055 98.39766 

MASK ANGLE- 1 DEG 


I 

EL EV 

(SINGLE, 

DUAL, 

CLEAR A) NLAT ELON 

TIME 

RANG 

1.000 

1.079 

3.592 

2.756 

2.219 -69.334 

263.313 

1.150 

2.989 

2.000 

5.205 

0.897 

0.702 

0.612 -72.583 

258.031 

1.167 

2.582 

3.000 

10.174 

0.440 

0.350 

0.315 -75.652 

250.604 

1.183 

2.179 

4.000 

16.527 

0.265 

0.213 

0.195 -78.396 

239.637 

1.200 

1.785 

5.000 

25.313 

0.172 

0.140 

0.130 -80.528 

223.122 

1.217 

1.408 

6.000 

38.572 

0.115 

0.095 

0.089 -81.569 

200.274 

1.233 

1.071 

7.000 

58.628 

0.081 

0.068 

0.065 -81.124 

175.823 

1.250 

0.825 

8.000 

68.386 

0.073 

0.062 

0.060 -79.380 

156.536 

1.267 

0.765 

9.000 

47.904 

0.095 

0.079 

0.075 -76.841 

143.511 

1.283 

0.929 

10.000 

31.266 

0.140 

0.114 

0.107 -73.880 

134.820 

1.300 

1.230 

11.000 

20.539 

0.212 

0.172 

0.158 -70.688 

128.769 

1.317 

1.591 

12.000 

13.127 

0.337 

0.269 

0.245 -67.364 

124.337 

1.333 

1.978 

13.000 

7.543 

0.604 

0.476 

0.423 -63.957 

120.938 

1.350 

2.378 

14.000 

3.036 

1.598 

1.238 

1.049 -60.496 

118.227 

1.367 

2.783 

15.000 

2.627 

1.868 

1.443 

1.212 -66.368 

242.138 

2.783 

2.835 

16.000 

6.994 

0.654 

0.515 

0.456 -69.714 

238.066 

2.800 

2.432 

17.000 

12.329 

0.360 

0.287 

0.260 -72.947 

232.587 

2.817 

2.036 

18.000 

19.251 

0.226 

0.183 

0.169 -75.988 

224.838 

2.833 

1.653 

19.000 

28.865 

0.151 

0.123 

0.115 -78.680 

213.341 

2.850 

1.298 

20.000 

42.619 

0.105 

0.087 

0.082 -80.714 

196.072 

2.867 

1.004 

21.000 

57.172 

0.082 

0.069 

0.066 -81.596 

172.672 

2.883 

0.837 

22.000 

52.827 

0.088 

0.073 

0.070 -80.981 

148.547 

2.900 

0.875 

23.000 

37.117 

0.119 

0.098 

0.092 -79.119 

130.029 

2.917 

1.097 

24.000 

24.921 

0.174 

0.142 

0.132 -76.517 

117.610 

2.933 

1.419 

25.000 

16.414 

0.267 

0.215 

0.197 -73.522 

109.294 

2.950 

1.786 

26.000 

10.144 

0.441 

0.351 

0.315 -70.311 

103.472 

2.967 

2.175 

27.000 

5.200 

0.898 

0.702 

0.613 -66.975 

99.182 

2.983 

2.574 

28.000 

1.081 

3.592 

2.756 

2.219 -63.561 

95.876 

3.000 

2.979 
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Elevation Angles at McMurdo 


The time variation of elevation of EOS or Landsat over McMurdo Sound may be seen with some elevation and 
attenuation plots. The first satellite pass over McMurdo may be seen in the lower left comer of Figure 3-1. 
Elevation is shown with data points, and the corresponding atmospheric attenuation of a hypothetical 27 GHz link 
is shown with a solid line. The data points are chosen at one minute intervals, and most passes are seen to be 
limited to less than 10 minutes. The long period between passes is taken out of the figure so that almost two days 
of satellite passes may be fitted onto the figure. The passes are seen to group as several high passes followed by a 
group of very low passes barely above the assumed 5 degree mask. The simulated elevation after 16 days is shown 
in the following figure (Figure 3-2). The figure is strikingly similar to that for the initial passes, but other methods 
should be used to confirm a true 16-day periodicity. 

Figure 3-3 shows the result of an autocorrelation of the entire 21.75-day elevation simulation with itself. The 
autocorrelation peaks at 1.0 when the series is cross correlated with itself for zero phase delay. It is seen to have a 
cyclic correlation with later passes. As the lag times get longer, the cross correlation is rarely over 0.4 until the lag 
time is set at 16 days (Figure 3-4) when the autocorrelation peaks at over 0.7. Longer phase differences also match 
the initial autocorrelation function (labeled Rx) well, but not perfectly. 

Elevation Probability Density Functions 

Elevation angle probability density functions (pdfs) have been interesting and important since the early ‘60s. Lutz 
and Dorosheski (9) recognized that satellite ground station economics depended on the ground station elevation 
angles. They looked at the entire ensemble of ground stations seen by a geosynchronous satellite, and treated the 
ensemble of elevation angles as a pdf. Sugai and Christopher (10) derived a simpler pdf which was useful for 
medium altitude to geosynchronous orbits. The simpler pdf was of the form: 


p(E) = 


f- 

1 180 j 

| COS 

— E 
180 , 

)('■ 

-2s sin 

f— £ i) 

U80 J) 


f 

1- 

< 

s 

-] 

36) 

-sin 
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Where 

E = elevation, deg. 

Emm = minimum elevation, deg. = 5° 

s = earth radius/semimajor axis = 1/6.6227 for geo satellite 


Figure 3-5 shows the emphasis on low elevation angle for the pdf generated by equation 3-1 for Geo altitudes. 
However, the 16-day simulation of Landsat4 passes over McMurdo Sound and shows an entirely different emphasis 
on low elevation angles, as seen in Figure 3-6. A fit for the McMurdo elevation pdf may be found as: 


p(E) = 0.00064269 csA— e] 
{ 180 ) 


0.003228 n cos 
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180 J 
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VISIBILITY TIME, MIN 


— Rx 

Figure 3-3: Elevation Autocorrelation vs. Lag Time, Min. 
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Figure 3-4: Elevation Autocorrelation Compared to 16-Day Autocorrelation 
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elevation, deg. 

Figure 3-5: Ground Elevation Probability Density Function to Geosynchronous Satellite 



0 20 40 60 80 

elevation, deg. 


Figure 3-6: Elevation PDF for McMurdo Sound-EOS 


Equation 3-2 is compared with the McMurdo data in Figure 3-7. It is useful, because it emphasizes the importance 
of the ground site’s ability to acquire Landsat or EOS quickly at the horizon. The overall communication time 
available increases sharply if the minimum elevation at the ground site can be 3 deg., for example, rather than 
5 deg. 


341 


WR96039.DOC 




Appendix A shows the full 16-day periodicity of actual Landsat4 elements (assumed to be identical to 
EOS elements here). The interarrival interval statistics and other link statistics may be summarized as in 
Table 3-2. 


Table 3-2: EOS Downlink Statistics to McMurdo Sound 
F = 8.28 GHz Precipitation Height = 1 .0 km 0.997 Availability 



Mean 

Standard Deviation 

Maximum 

in view, hr (min) 

0.1848(11.08) 

0.032 (1.92) 

0.2167 (13.0) 

out of view 

1.4609(87.6) 

0.0526(3.15) 

1.550(93.0) 

attenuation, dB 

0.70 

0.80 

3.59 

elevation (deg) 

17.653 

16.868 

90 

Rate (deg/sec) 
elevation 

0 

0.1353 

0.586 

azimuth 

-0.057 

0.3495 

18.82 


With maximum antenna rates from a separate run with 5 sec. increments 


4.0 Outages For General Sun Synchronous Satellites 

A large range of altitudes of sun synchronous satellites may be considered for use with a McMurdo Sound ground 
site. Low-altitude, low-cost satellites such as reported for a November 1994 Pegasus launch may be considered a 
prime opportunity for the McMurdo site. This section shows the variation in outage statistics as a function of 
altitude. 

The inclination of sun synchronous satellites is actually function of altitude (Figure 4-1). This feature, along with 
the increased viewing region associated with high altitudes, has noticeable effects on outage times and viewing 
times. Long simulation can be done at each altitude and outage statistics extracted. Figure 4-2 shows the mean 
outage times ranging from 1.45 hours at 550 km altitude to 1.5 hours at altitudes slightly greater than 1000 km. 
Altitudes much less than 500 km, however, may actually skip an orbit before seeing McMurdo again. This appears 
as long outage time in Figures 4-2 and 4-3. 
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Figure 4-1: Required Inclination for Sun Synchronous Satellites 



Figure 4-2: Outage PDF vs. AititudelOOO 
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0.4 0.6 0.8 1 1.2 

Altitude (+ 1000) km 

Figure 4-3: Contour Plot of Probability Density Function 


Conclusions 

We have discussed outage intervals, innage intervals, and dynamic propagation losses for a sun synchronous- 
McMurdo Sound ground site. Antenna rates have also been extracted from the link dynamics. McMurdo Sound 
has been shown to be an attractive ground site for a wide range of sun synchronous satellites from 500-1000 km 
altitude. EOS and Landsat were examined in detail. The observed 16-day Landsat periodicity was closely 
reproduced in the simulation. A unique elevation probability density function for McMurdo Sound was shown. 
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APPENDIX A: CONDENSED ELEVATION ANGLE 
RECORD AT MCMURDO SOUND 


A tabular output for ground antenna angles, times, and propagation loss is indispensable, and can be obtained from 
an April 1995 STel report. The figures below extract all the outage time to allow almost 1.7 days of EOS passes to 
be presented on each figure. They also allow a hypothetical 27 GHz downlink to be examined for atmospheric losses 
(1 1). Figure A-l shows elevation angle vs time for the first 1.7 days. A quasi-cyclic period can be seen every day 
when very low elevation angles are seen at McMurdo. The elevation angles (data points) are read off the. left scale. 
Atmospheric attenuation is read off the right scale. 

Figures A-2 to A- 10 extend the elevation angle record consecutively to 17.28 days. 


ELEVATION FOR McMURDO-EOS 
FOR SUCCESSIVE PASSES 


McMUROO PASSEStOAYS 1.7- 3.3 


ELEV, DEG A ciB 



0 50 100 150 200 250 

TIME, MIN 

EIEVATION.DEG AT TEN ,<JR 


EXAMINE t7 QHx LINK LOSS AT P«0,M7 



240 260 200 300 320 340 360 380 400 4 20 4 40 460 4fln 
T.MIN FOR EACH PASS 

— “ ELEV AT TEN .dB 


A-1 


A-2 


MoMUROO PR8SE84MY8 3.S-6.1 

DEG dB 



400 500 520 540 560 580 600 620 640 660 680 700 7?n 
T.MIN FOR EACH PASS 

ELEV ATTFN.dR 


MqMURDO PASSE 3;0 AYS 6.1-6.90 



720 740 760 700 800 820 040 860 880 900 920 940 900 
T.MIN FOR EACH PASS 

^ ELEV AT TEN. dB 
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A-4 
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McMUROO PA88ES;DAY8 6.08-8.63 


McMURDO PASSE8;DAY3 6.63-10.34 


0 96 0 98 

1 1 02 1 04 1 06 1 08 11 1 12 114 116 1 18 1 ? 

T.MIN FOR EACH PASS (Thousands) 

1.2 122 124 126 128 13 1 32 1 34 136 138 14 

T.MIN FOR EACH PASS (Thousands) 


— — ELEV ATTEN .(IB 

~ELEV ATTEN.dB 


4-5 

A-6 


McMURDO PASSES;DAYS 10.34-12.07 

McMURDO PASSES;DAYS 12.1-13.9 DAYS 

DEG 

nn 

dP 

DEG 

80 1 


144 146 148 1.5 1.52 *154 156 158 16 16? 164 166 160 

T.MIN FOR EACH PASS (Thousands) 

168 17 172 1.74 176 1 78 18 1 82 1 84 1 86 1 88 19 

T.MIN FOR EACH PASS (Thousands) 

— — ELEV ATTEN.dB 

— ELEV ATTEII.rtO 

A-7 

4-0 

McMUROO MMESsDAYB 13.9-1S.M DAYS 

McMURDO F*88E8|DAY8 16.66-17.28 

DEG 0B 

DEG. 

mn — 


192 194 1.96 198 2 2 02 2 04 2 06 2 08 2 1 2 12 2 14 2 16 

T.MIN FOR EACH PASS ( Thousands) 


2 16 2 18 2 2 2 22 2 24 2 26 2 20 2 3 2 3? 2 34 2 36 

T.MIN FOR EACH PASS (Thousands) 


‘ ELEV ATI FI I OH 


A-10 
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